Low energy collective excitations in a superfluid trapped Fermi gas 
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ON ' We study low energy collective excitations in a trapped superfluid Fermi gas, that describe slow 

variations of the phase of the superfluid order parameter. Well below the critical temperature the 
corresponding eigenfrequencies turn out to be of order the trap frequency, and these modes manifest 
themselves as the eigenmodes of the density fluctuations of the gas sample. The latter could provide 
an experimental evidence of the presence of the superfluid phase. 
PACS number(s): 03.75.Fi, 05.30.Fk 

The search of a novel behavior of a gaseous system in the confined geometry of a trap attracts a lot of attention since 
the discovery of BEC in trapped gas samples of alkali atoms (l]-[| . This behavior is mostly related to the combined 
• effects of an interparticle interaction and quantum statistics, which plays a dominated role below the degeneracy 
\ temperature. In the case of a Bose system, the latter corresponds to the critical temperature of BEC. For a Fermi 
system the degeneracy temperature equals to the Fermi energy Tp = sp, and at temperatures below Tp the Pauli 
exclusion principle should be taken into account . Remarkably, that at much lower temperature (called the critical 
temperature T c -C Tp) there could be a transition into a superfluid phase, as a result of Cooper pairing. Possible 
versions of the transition have recently been discussed in Refs. Although the pairing strongly influences only 



^ . a small fraction (~ T c /ep <C 1) of quantum states in the vicinity of the chemical potential, it nevertheless changes 



C3 



the behavior of the Fermi system significantly, because precisely those states govern the response of the system to a 
small external perturbation. For example, as it was shown in Ref. JHJ], the superfluid pairing changes (smears out) 
the resonance structure of the gas density oscillations in a parabolic trap. 

In this paper we study low energy collective excitations in a superfluid Fermi gas trapped in an isotropic harmonic 
potential. These excitations correspond to the phase fluctuations of the order parameter A(r) and, as shown, well- 
below the transition temperature their eigenfrequencies are of the order of the trap frequency. The damping of the 
excitations is expected to be small, and, therefore, these modes are well-defined collective excitations. They manifest 
themselves as oscillations of the superfluid component of the Fermi gas, and hence, at low enough temperature, when 
, the density of a normal component tends to zero exponentially, these collective excitations define the eigenmodes for 
^ ■ the density fluctuations of the entire gas sample. As these modes exist only below the transition temperature, their 
OO experimental observation could provide an indication of the presence of the superfluid phase. 

Specific properties of the collective modes in a superfluid phase depend on the structure of the order parameter and, 
hence, on the type of pairing. For a singlet s-wave pairing the order parameter is a complex function, and one has 
two branches of collective excitations corresponding to the phase and modulus variations of the order parameter. For 
a triplet p-wave pairing the order parameter is a complex 3x3 matrix, and, hence, additional branches of collective 
Q\ [ modes appear (see, e.g., Ref. However, in both cases the lowest in energy branch corresponds to the fluctuations 

■ of the phase of the order parameter (Bogolyubov sound). In this paper we study this mode for a trapped superfluid 
Fermi gas. For simplicity we consider the case with a "singlet" s-wave pairing, that implies the presence of an 
attractive interatomic interaction. This situation can be realized in a gas of 6 Li atoms in a magnetic trap, where 
they are characterized by a large and negative triplet s-wave scattering amplitude a » — 1 140 A jl4j between the two 
hyperfine components. 

We consider a two component gas of fermionic atoms (a and [3 atoms) trapped in an isotropic harmonic potential. 
We assume the two components (for example, hyperfine components) have equal masses and concentrations. The 
Hamiltonian of the system is (h = 1) 



-a 
c 
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H = E / dv4{v)H^{v) + V J dr^(r)^ Q (r)^(r)^(r), 



(1) 
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where tpi(r) with i = a, (3 are the field operators of a and (3 atoms, Hq = —V 2 /2m, + mfl 2 r 2 /2 — fi, CI the trap 
frequency, m the mass of the atoms, and // the chemical potential. The second term in Eq. ([j]) corresponds to an 
attractive elastic short-range interaction between a and f3 atoms (s-wave scattering length a < 0) with V = Aita/m 
being the coupling constant. In the Hamiltonian ([!]) we neglect the interaction between atoms of the same sort, 
originating in the case of fermions only from the scattering with orbital angular momentum I > 1. The presence of 
an attractive intercomponent interaction in the s-wave scattering channel leads to a superfluid phase transition j^] 
with the critical temperature T c <C /i. We assume that T c is much larger than fl and, hence, the value of the critical 
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temperature in the trap is very close jlO| to the critical temperature in a spatially homogeneous gas with the 
density being equal to the maximum density n of the trapped gas sample, = 0.28£i?exp{— 1/A} jl5|], where 
A = 2\a\p F /it is the small parameter of the theory (A < 1), pp = (37r 2 n ) 1 / 3 and £p = p 2 F /2m. 

The superfluid phase is characterized by the order parameter (complex function) A(r) = \V\ (tp a (r)jp^(r)) . The 
equilibrium form of A(r) for a trapped gas sample was studied analytically in Ref. O] and numerically in Ref. [ p7[ . 
The appearance of A(r) strongly influences only the quantum states within the range of order T c near the chemical 
potential level. As a result, the gas density profile changes only slightly (An(r)/n(r) ~ T c /n -C 1), and hence, 
one has n(r) = rio(l — (t/Rtf) 2 ) 3 ^ 2 m the Thomas-Fermi approximation for both normal and superfluid phases 
(Rtf — PF/mQ is the Thomas- Fermi radius of the gas cloud). The interparticle interaction leads to corrections to 
this formula. However, the leading the mean-field corrections ]l6| ] are proportional to the small parameter of the 
theory A, and hence, will be neglected. 

At finite temperatures it is convenient to study the evolution of the system in imaginary (Matsubara) time r € 
[0, 1/T] jL8). If one splits the order parameter A(r,r) into its equilibrium part Ao(r) (Aq(i-) = Ao(r)) and a small 
fluctuation S(r, r) = T ^ w <5 w (r) exp(— iwr), where u> = 7rT(2n + 1) with n being an integer, is the Matsubara 
frequency, then the equation for <L(r) reads (see, e.g. [?] ) 



\v\ 



= TJ2J dv'{G^+^y)G-^y)^{v') - F^iv^F^v^Sl (r')} (2) 



(the field 8*(r) obeys a complex conjugate equation). In Eq. (g) G w (r,r'), F^(r, r') are the normal and anomalous 
Green functions of the Hamiltonian 



i—ct,f3 

in the Matsubara representation. They can be expressed through solutions (u v (r),v^(r)) of the Bogolyubov-de Gennes 
equations 

:S) +a » w (-"$>)■ < 3 » 

which define the Bogolyubov (unitary) transformation from the bare fermions tp a , ipp to single-particle excitations 
a v , /3„with energies e u > 



Mr) \ _ 
#( r ) 



E 



(4) 



such that the Hamiltonian Hq is diagonal. Then the Green functions are 

G.( ri ,r 2 ) = Y ( M - (riK(r2) + ^ (riK(r2) l 



FUn,*) = £ l uAriK{r2) - <(riK(r2) | . (5) 

' I IUJ - E v IUJ + E v I 

For the fluctuation of the phase of the order parameter one has 5 w (r) = — <5* (r) = 2iAo(r)<p [1J (r) with </J w (r) <C 1 
being a real function, and Eq. (0) can be rewritten in the form 



dr'^ (r,r')^(r') = dr'K w (r,r')5 u (r'), (6) 



where X = T^ Wl (G W1 G_ W1 + F U1 F U1 ) is the static kernel and X w = TY, Ul { G ^+^ G -^i + ~ the 

dynamic one. 

For (T c — T)/T c <C 1 the order parameter is small, A(r) <C T c , and one can neglect the contribution from the 
anomalous Green functions. Under the assumption T c ^S> 0, the normal phase Green functions G w (r, r') can be taken 
in the semiclassical approximation. Then for a slow varying function ip one gets the following results for the left and 
right hand sides of Eq. (||) 



2 



W\ 



r- / dr'K (r,r')5 u (r') = 2iN(R) 
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(1 _ ij2)3/2 Vay , + 2(l-i^)V R A V R ^) , (7) 



where £(2) is the Riemann zeta function, N(R) — (mpp /ir 2 )\/l — R 2 the density of states on a local Fermi surface 
(R = v/Rtf), and 



J dv'K^v,v')5^v') = 2iA Q (R)N(R) 
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For T <§; T c the calculations are more laborious. In this case one can use, for example, Eqs. (g) and various relations 
between u v and v v , that follows from Eqs.. (Q) and from unitarity of the Bogolyubov transformation, Eq. (^). The 
answer for the left hand side of Eq. (o) is 
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The kernel K^r^r') can be calculated by usin g E qs. (j^) with the semiclassical solutions of the Bogolyubov-de 
Gennes equations (|3|), which were found in Ref. |1 1| for the considered case of a spherically symmetric harmonic 
trapping potential. In this case there are two different types of excitations: in-gap and above-gap. The eigenenergies 
of the in-gap excitations are smaller than the maximum value of the order parameter A(0), and their wave functions 
(w„(r), v v (r)) are pushed by A(r) to the outer part of the gas sample Be. For this reason these states give only 
exponentially small contribution (~ exp(— T c /fi)) to the right hand side of Eq. (^), and, hence, only above-gap 
excitations are important. The calculations with slow varying function ip then yield 



y"dr'i^(r,r')<Ur') = -2i 



2 N(R) 



(10) 



(This result can also be obtained by using the Green functions G w (r, r') and F w (r, r') in the local density approxima- 
tion.) 

After making an analytic continuation in Eqs. (@) and ([l0]) from the Matsubara frequency w to a real one E, and 
combining the results with Eqs. ^) and (|]), we finally obtain the equations describing the collective modes related 
to the phase fluctuations of the order parameter (Bogolyubov sound). 

For (T c - T) /T c < 1 one has 



7Q 2 C(3) 
6vr 3 T c \VT^W 



R 



(1 - i? 2 ) 3/2 V R (p(R)l + 2(1 - i? 2 )V R ln AoVr^(R) ) = iEip(K) 



(11) 



and it follows from this equation that for temperatures close to the critical temperature T c the eigenenergies E are 
purely imaginary: the collective mode decays rapidly into pairs of single-particle excitations. 
For T <^T C the equation reads 
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(1 _ Jj2)3/2 VR p( R ) =E 2 l f(R), 
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and in this case the eigenenergies E of the collective modes are real and of order the trap frequency f2. These modes 
being excited results in oscillations of the superfluid current j = (n/m)V<p and density n = uq + Sn, that are related 
to each other by the continuity equation dn/dt + divj = 0. (For T -C T c one can neglect the contribution of the 
normal component of the gas and hence, j s = j and n s — n.) As a result, the entire gas sample oscillates according to 
the formula 



n(r, t) — n (r) + <5n(r, t) 



1 + — vV 

m 



1 



n (r -\ VV>) 



where ip(r,t) — j 1 <p(r, t')dt'. (The prefactor in the square brackets ensures the conservation of a total number of 
particles.) 

The damping of the collective modes, which does not contained in Eq. ( |l2| ) , will be mainly determined by processes 
of decay and scattering on in-gap excitations (these mechanisms of damping of collective modes in Bose-condensed 
atomic gases are often called as Beliaev and Szepfalusy-Kondor mechanisms respectively; see, e.g. pf|). One can say 
that the energy of the collective mode converts into a normal component generated in the outer part of the gas sample. 
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But because the coupling between the order parameter fluctuations and the in-gap excitations is exponentially weak 
(~ exp(— T c /f2)) , the damping is expected to be small. 

It should be mentioned that Eq. ( fl2"| ) can be obtained in the framework of a hydrodynamic description of a 
superfluid Fermi gas. When the superfluid velocity v s = m V<£ and the deviation of the particle density Sn from 
its equilibrium value no are small, the corresponding Hamiltonian has the form 



2 m 



U (no + Sri) 



2 m 
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U"{no)Sn 2 + [/(n ) 



(13) 



where U(ri) is the density dependent part of the energy. (The equilibrium density no is defined by the condition 
{/'(no) = 0.) In the Thomas -Fermi approximation, 



U(ri) 



p 2 dp 

P < PF (n) 2m (2tt) 3 
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H n, 



(14) 



and the equilibrium gas density profile no(r) — (p|,/37r 2 )\/l — R 2 , as it should be. (In Eq. ( |14| ) we do not include the 
effects of the mean field interaction and the superfluid pairing because they contain small parameters A and (T c /ef) 2 
respectively.) For the quantity U"(n ) in Eq. ( |l3| ) one now has {/"(no) = (3ir 2 )~ 2 / 3 N(r)~ 1 , and the standard 
commutation relation [6n(ri), <^(r 2 )] = iS(ri — r 2 ) leads to the equations 

dip/dt = i[H h ,ip}= U"(n )Sn 
d(Sn)/dt = i [H h , Sn] = -V(n V^)/m, 



(15) 



from which one immediately gets Eq. (pL 2|) for the phase dynamics and the equation 



d 2 (Sn)/dt 2 + ^-V R 



(1 - i? 2 ) 3 / 2 V R 



Sn 



= o 



for the density fluctuations. 

Equation (jlj) (or (E^)) together with the condition that ip (or Sri) is finite, provides us with the eigenfrequencies 
and eigenfunctions of the collective modes of the superfluid Fermi gas. For a spherically symmetric breathing modes 
(angular momentum I is zero) one has for the eigenfrequencies 



{E n o/ttf 



-n(n + 2), n = 1,2,.. 

o 



(ri = corresponds to a constant phase and, hence, has no physical significance), and for the eigenfunctions 

3 

</? n o(R)oc 2 F 1 (-n,n + 2;-;R 2 ), 

where 2^1 (« 7 b; c; z) is the hypergeometric function. 
For nonzero angular momentum 



(16) 



(17) 



{E nQ /n) 2 = l + -n(n + 1 + 2), n = 0, 1, 2, . 



(18) 



and 



<p nl (R) cx R l 2 F 1 (-n, n + l + 2-- + l- R 2 )Y lm {6, , 



(19) 



(Note, that the eigenfunctions (|l7j ) and ( |l9|) are orthogonal with the weight l/\/l — R 2 .) 

The spectrum (16) and ( p^ ) coincide with that for a trapped normal Fermi gas in a hydrodynamic regime 
and the lowest eigenfrequencies are equal to those for a classical gas |2^]. The lowest eigenmodes are of a special 
interest because they can be excited by a modulation of the trap frequency. (A small external perturbation V ext 
results in an extra term — iEV e *t: or dV cyi t/dt in the time representation, on the right hand side of Eq. (UJ).) It is 
reasonable to assume that the trapped Fermi gas just above T c is in a collisionless regime (fir ~ Q(\ 2 T^Je f) _1 ~ 
(il/X 2 T c ) exp(l/A) 3> 1), where the lowest monopole, dipole and quadrupole collective modes were calculated in 
p3[ . For the superfluid phase, as follows from Eqs. ( |l6| ) and (p^), the lowest eigenfrequency E\o for the monopole 
breathing mode ({ = 0, n = 1) is equal to 2fl (this result can be obtained on the basis of the sum rules p3), and 
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one has the anticipated result Eqi = f2 for a dipole mode I = 1, n = 0. These eigenfrequencies coincide (within 
the considered accuracy) with those from Ref. p3[ . On the other hand, Eq. ( p"8| ) gives E02 — \Z2fl for the lowest 
quadrupole mode, while one has E02 — 2fl [^3| in the collisionless regime above T c . Experimentally the quadrupole 
mode can be excited by a small anti-phase modulation of the trap frequency in, for example, x and y directions, 
K x t(r, t) = (mfl 2 /2)(x 2 — y 2 )( cos(wi) with £ -C 1, and the response of the gas sample for T > T c and T<T C will 
have resonances at different frequencies, 217 and ^/2H, respectively. 

In conclusion, we have found the low energy collective modes in the superfluid trapped Fermi gas. These modes 
are related to the fluctuations of the phase of the superfluid order parameter, and, hence, describe the motion of 
the superfluid component of the gas. Just below the critical temperature of the superfluid phase transition the 
eigenenergies of the collective modes are purely imaginary, and these modes describe the diffusive relaxation of a 
superfluid fluctuation. For temperatures well below T c , the eigenenergies are of order the trap frequency, and the 
damping is small. Therefore, these modes can manifest themselves as the eigenmodes of the gas density oscillations 
that can be observed experimentally, and serve as an indication of the superfluid phase transition. 
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M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, and E.A. Cornell, Science, 269, 198 (1995). 

C. C. Bradley, C.A. Sackett, J.J. Tolett, and R.G. Hulet, Phys. Rev. Lett., 75, 1687 (1995). 

K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Durfee, D.M. Kurn, and W. Ketterle, Phys. Rev. Lett., 
75, 3969 (1995). 

Th. Busch, J.R. Anglin, J.I. Cirac, and P. Zoller, e-print |cond-mat/9805037 . 
L. Vichi, M. Inguscio, S. Stringari, G.M. Tino, J. Phys. B 31, L899 (1998). 
B. DeMarco and D.S. Jin, Phys. Rev. A 58, R4267 (1998). 

H.T.C. Stoof, M. Houbiers, C.A. Sackett, and R.G. Hulet, Phys. Rev. Lett. 76, 10 (1996). 
MA. Baranov, Yu. Kagan, and M.Yu. Kagan, JETP Lett. 64, 301 (1996). 

M. Houbiers, R. Ferwerda, H.T.C. Stoof, W.I. McAlexander, C.A. Sackett, and R.G. Hulet, Phys. Rev. A 56, 4864 (1997). 
MA. Baranov and D.S. Petrov, Phys. Rev A 58, R801 (1998). 
MA. Baranov, preprint |sond-mat/9801142| . 

D. Pines and Ph. Nozieres, The theory of quantum Liquids, (W.A. Benjamin, INC. New- York Amsterdam, 1966). 

D. Vollhardt, P. Wolfle, The superfluid phases of helium 3, (Taylor&Francis, London, 1990). 

E. R.I. Abraham, W.I. McAlexander, J.M. Gerton, R.G. Hulet, R. Cote, and A. Dalgarno, Phys. Rev. A 55, R3299 (1997). 
L.P. Gor'kov and T.K. Melik-Barkhudarov, Zh. Eksp. Teor. Fiz. 40, 1452 (1961) [Sov. Phys. JETP 13, 1018 (1961)]. 
G.M. Bruun and K. Burnett, Phys. Rev. A 58, 2427 (1998). 

G.M. Bruun, Y. Castin, R. Dum, and K. Burnett, Euro. J. Phys. (in press). 

L.D. Landau and E.M. Lifshitz, Statistical Physics, Part 2, (Pergamon Press, Oxford, 1980). 

G. Eilenberger, Zeitschr. fur Phys. B 214, 195 (1968). 

P.O. Fedichev, G.V. Shlyapni kov, and J.T.M. Wa lraven, Phys. Rev. Lett. 80, 2269 (1998). 
G.M. Bruun and C.W. Clark, |eond-mat/9905263| . 

D. Guery-Odelin, F. Zam belli, J. Dalibard, S . Stringari, |cond-mat/9904409 
L. Vichi and S. Stringari, cond-mat/9905154. 
A.J. Leggett, private communication. 



5 



